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Abstract 

To determine scale velocity an additional similarity 
parameter is needed to supplement the Ruff scaling 
method. A Weber number based on water droplet MVD 
has been included in several studies because the effect 
of droplet splashing on ice accretion was believed to be 
important, particularly for SLD conditions. In the pre- 
sent study, ice shapes recorded at Appendix-C condi- 
tions and recent results at SLD conditions are reviewed 
to show that droplet diameter cannot be important to 
main ice shape, and for low airspeeds splashing does 
not appear to affect SLD ice shapes. Evidence is pre- 
sented to show that while a supplementary similarity 
parameter probably has the form of a Weber number, it 
must be based on a length proportional to model size 
rather than MVD. Scaling comparisons were made be- 
tween SLD reference conditions and Appendix-C scale 
conditions using this Weber number. Scale-to- 
reference model size ratios were 1:1.7 and 1:3.4. The 
reference tests used a 91-cm-chord NACA 0012 model 
with a velocity of approximately 50 m/s and an MVD of 
160pm. Freezing fractions of 0.3, 0.4 and 0.5 were 
included in the study. 

Nomenclature 

A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

c Airfoil chord, cm 

c p Specific heat of air, cal/g K 

c P ,ws Specific heat of water at the surface tempera- 
ture, cal/g K 

d Cylinder diameter or twice the leading-edge 

radius of airfoil, cm 

h c Convective heat-transfer coefficient, 

cal/sm 2 K 

h G Gas-phase mass-transfer coefficient, g/sm 2 
K Inertia parameter, dimensionless 

K 0 Modified inertia parameter, dimensionless 

L Length proportional to model size, cm 


LWC Cloud liquid-water content, g/m 3 
MVD Water droplet median volume diameter, pm 
n Freezing fraction, dimensionless 

P General similarity parameter, dimensionless 

p Pressure, Nt/m 2 

p w Vapor pressure of water in atmosphere, Nt/m 2 

p ww Vapor pressure of water at the icing surface, 
Nt/m 2 

r Recovery factor, dimensionless 

Re Reynolds number of model, dimensionless 

Re s Reynolds number of water droplet, dimen- 

sionless 

SLD Super-cooled large droplet 

tf Freezing temperature, °C 

t s Surface temperature, °C 

t Temperature, °C 

T Absolute temperature, K 

V Air velocity, m/s 

We Weber number based on droplet size and water 
properties, dimensionless 

We c Weber number based on model size and air 
properties, dimensionless 

We L Weber number based on model size and water 
properties, dimensionless 

po Collection efficiency at stagnation line, dimen- 

sionless 

(p Droplet energy transfer parameter, °C 

A Droplet range, m 

A-stokes Droplet range if Stokes Law applies, m 

Af Latent heat of freezing, cal/g 

A v Latent heat of condensation, cal/g 

p Air viscosity, g/m s 

6 Air energy transfer parameter, °C 

p Air density, g/m 3 

Pi Ice density, g/m 3 

p w Liquid water density, g/m 3 

a Surface tension of water over air, dyne/cm 
r Accretion time, min 

Subscripts 

R Reference 

5 Scale 

st static 

tot total 
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Introduction 
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VMVDp 
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(3) 


Various icing scaling studies over the past 50 years 
(see, for example, Ruff 1 ) have accepted that the most 
important similarity parameters affecting ice shape are 
the freezing fraction, n, accumulation parameter, A c , 
and modified inertia parameter, K 0 . The leading-edge 
collection efficiency, p 0 , is directly dependent on K n . 
For scale tests, these three parameters are matched to 
their respective reference values (established from the 
hill-scale or reference conditions to be simulated). A 
fourth parameter, either the droplet energy transfer pa- 
rameter, (p. or air energy transfer parameter, 8, was in- 
cluded in Ruff s scaling analysis. The four equations 
that result can then be solved for the scale temperature, 
LWC, time and MVD. For glaze ice accretions, an addi- 
tional parameter is needed to determine scale velocity. 
The identification of this parameter has been the subject 
of ongoing studies for several years. 2 3,4 Flere evidence 
from past studies along with some recent SLD ice 
shapes will be examined to identify the probable form 
of this parameter. Results from using this parameter in 
SLD-to-Appendix C scaling tests will be presented. 


Scaling Equations 

The basic defining equations for the similarity parame- 
ters used here followed Ruff. 1 The scaling method in- 
volved matching scale and reference values of the simi- 
larity parameters, A c and n, with other parameters 
evaluated for comparison. In Ruff s method, the scale 
velocity was set at the user’s convenience; in this study 
it was determined by matching a Weber number to be 
defined later. The equations for the remaining similar- 
ity parameters will be presented here. 

The modified inertia parameter, K 0 , was defined by 
Langmuir and Blodgettf 



A 


A,, 



(1) 


In equation (1), K is the inertia parameter, 

K Pw MVD 2 V 

18 dp 


(2) 


where d is the diameter for cylindrical models or twice 
the leading-edge radius for airfoils. For the NACA 
0012, a leading-edge radius of 0.0158c was used (see 
Abbott and von Doenhofr), where c is the model 
chord. AJAstokes is the droplet range parameter, defined 
as the ratio of actual droplet range to that if Stokes drag 
law for solid-spheres applied. It is a function only of 
the droplet Reynolds number, Res- 


Langmuir and Blodgett’s tabulation of the range pa- 
rameter was fit to the following expression for this 
study: 


1 

2 Stokes 


'0.8388 + 0.001483 Re^V 1 

v+ 0.1847y^ 


(4) 


Of more practical interest than K„ is the collection effi- 
ciency at the leading edge, p 0 , which was shown by 
Langmuir and Blodgett to be a function only of K 0 , 

,.84 
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The accumulation parameter is: 

LWCVt 


(5) 


( 6 ) 


If all the water impinging on the leading edge freezes at 
that location and the leading-edge collection efficiency 
is 100%, A c is a measure of the thickness of ice that will 
accrete relative to airfoil size. 


The freezing fraction is the ratio of the amount of water 
entering a specified region on the surface that freezes 
there. From Messinger’s 7 surface energy balance, the 
freezing fraction is 


n = 


A 




(7) 


The individual terms in this expression are (j>. the water 
energy transfer parameter, 

V 2 

<P = tf-t st -— ( 8 ) 

p,ws 


8, the air energy transfer parameter, 
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and b, the relative heat factor, introduced by Tribus, et. 
al. 8 


LWC V c p 
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(10) 
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Equation (9) from Ruff includes compressibility effects. 
A simpler form without compressibility was used by 
Charpin and Fasso 9 and others. Ruff s expression for 0 
was used in this paper, but values found without com- 
pressibility are not significantly different for most icing 
conditions. 

Rationale for Development of Additional Parameter 

Because the Ruff method does not restrict the value of 
scale velocity, an additional similarity parameter can be 
used to determine V s . Evidence that such a parameter 
was needed was provided by Bilanin and Anderson" in 
studies of the effect of changing surface tension of wa- 
ter. Figure 1 includes published and unpublished data 
from their studies and compares ice shapes with and 
without surfactant addition to the spray for velocities of 
47, 67 and 94 m/s. The shaded ice shapes were ob- 
tained with the NASA Glenn Icing Research Tunnel 
(IRT) spray bar system de-mineralized water supply 
with no surfactant added, and the shape represented by 


the solid line resulted from surfactant addition to the 
water. Other than the addition of surfactant the same 
test conditions were used for each pair of tests, as 
shown in the accompanying table. Surfactant addition 
reduced the surface tension to roughly half that of wa- 
ter. Because n, and A c were the same for each pair of 
tests, the leading-edge ice thickness was also nearly the 
same. Flowever the included angle between the horns 
decreased dramatically when surfactant was added. 
Florn angle also decreased when velocity increased as 
can be seen by comparing figure 1(a) with (b) and (c). 
All the results shown were made at approximately the 
same freezing fraction. Thus, both surface tension and 
velocity have an effect on ice shape independent of 
freezing fraction, and scale velocity cannot be chosen 
arbitrarily. Clearly, then, a similarity parameter de- 
pendent on the ratio V“lo h must be included in scaling 
methodology, where the powers a and b are not yet 
determined. 

Additional evidence of the form of this parameter can 



JfC 


x c 


x/c 


Water, 6-9-94 Run 4 


Water. 8-27-93 Run 9 Water. 8-27-93 Run 8 


Water + Surfactant, 
8-30-93 Run 4 


Water + Surfactant. 
8-30-93 Run 5 


Water + Surfactant. 
8-30-93 Run 3 


(a) V= 47 m/s 





(b) 

V=67 

m/s 




(c) V 

= 94 

m/s 




Date/Run 

d, 

cm 

°c 

V, 

m/s 

MVD, 

pm 

LWC, 

g/m 3 

min 

a 

dyne/ 

cm 

fio, 

% 

A c 

n 

b 

<t>, 

°C 

e, 

°c 

Re, 

10 4 

We, 

10 3 

We c , 

10 3 

We h 

10 6 

(a) 6-9-94/4 

5.1 

-7.8 

46.9 

40.0 

1.16 

16.0 

65 

66.0 

1.13 

0.28 

0.76 

7.6 

11.3 

18.5 

1.36 

2.23 

1.73 

8-30-93/4 

5.1 

-8.0 

46.8 

40.0 

1.17 

16.0 

30 

66.0 

1.13 

0.29 

0.76 

7.7 

11.5 

18.5 

2.91 

4.80 

3.71 

(b) 8-27- 

5.1 

-11.8 

66.9 

34.1 

1.39 

10.2 

65 

65.4 

1.22 

0.33 

1.08 

11.3 

15.7 

26.7 

2.35 

4.53 

3.51 

8-30-93/5 

5.1 

-11.8 

67.1 

34.1 

1.39 

10.2 

30 

65.5 

1.23 

0.32 

1.08 

11.3 

15.6 

26.8 

5.12 

9.88 

7.66 

(c) 8-27-93/8 

5.1 

-12.2 

94.0 

30.0 

1.10 

9.0 

65 

65.6 

1.20 

0.31 

1.03 

11.2 

14.0 

36.6 

4.08 

8.7 

6.94 

8-30-93/3 

5.1 

-12.0 

93.6 

30.0 

1.10 

9.2 

30 

65.6 

1.22 

0.31 

1.03 

10.9 

13.7 

36.4 

8.75 

18.7 

14.88 


Figure 1 . Effect of Surfactant and Velocity on Ice Shape. Vertical Cylinders Tested in the NASA Glenn IRT. Pub- 
lished and Unpublished Ice-Shape Data from Tests by Bilanin and Anderson. 


NAS A/CR— 2005-2 13850 


3 



be found from tests performed in the IRT in 1998 by 
Chen. 10 Chen studied a 61-cm-chord GLC 305 airfoil 
model in which only the droplet MVD was varied. Re- 
sults for droplet sizes of 55 and 20pm at freezing frac- 
tions of 0.3 and 0.5 are given in figure 2. At each freez- 
ing fraction, the values of A„ n, b, (/>, 6 and Re were 
maintained constant as drop size decreased. Although 
Po changed from about 92 to 74%, there was no meas- 
urable effect of reducing droplet size on the main ice 
shape. Undoubtedly, the icing limit would have 
changed with droplet size (that is, with /%), but this 
feature was not measured. Chen’s results, and others 
like them for NACA 0012 airfoils, show that the pa- 
rameter being sought cannot be dependent on MVD. 
The present study included tests to evaluate if Chen’s 
observations were also true for SLD conditions, and the 
results will be presented below. 

Chen’s investigation also included an evaluation of the 
effect of model size. Figure 3 illustrates that with p n , 
A c , n, b, (f> and 6 constant a reduction in chord from 91 
to 30 cm moved the glaze horns rearward. This is the 
same effect shown in figure 1 for decreasing velocity. 

Previous studies with both cylinders and NACA 0012 
airfoils 11, have shown that the temperature and LWC 
do not have effects on the ice shape independent of the 


freezing fraction. Thus, the general form of the sup- 
plementary parameter must be 

V“c c 

P = const — — (11) 

a b 


This form suggests a Weber number based on chord: 


We c 


V 2 cp 

a 


( 12 ) 


Studies by Bartlett 1 ’’ 14 and Oleskiw, et. al. 15 found no 
measurable effect of pressure on ice shape. These ob- 
servations rule out the dependence on air density in 
equation (12), making water density a better choice. 
Furthermore, the length may not be chord itself but 
rather some physical characteristic related to the accret- 
ing ice that is proportional to chord; for example, the 
water-film thickness. Because this length is not yet 
identified, L will be used to represent it, and the Weber 
number is then 


We, 


V 2 Lp n . 

o 


(13) 


The trends apparent in figures 1 and 3 show that an 
increase in We L has the same effect on glaze ice shape 
as an increase in /?; that is, horns move back with an 
increase in the included horn angle as either n or We L 
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pm. Run 924.43 



(a) n = 0.28 








(b) 

n = 0. 

52 








Run 

c, 

cm 

Lh 

°c 

v, 

m/s 

MVD, 

pm 

LWC, 

g/m 3 

T, 

min 

Pn, 

% 

A c 

n 

b 

4>, 

°C 

e, 

°C 

Re, 

10 4 

We, 

10 3 

We 0 

10 3 

We L , 

10 6 

(a) 

624.40 

61.0 

-10 

90 

55 

1.16 

6.1 

91.6 

2.53 

0.28 

0.84 

8.9 

11.3 

10.9 

6.76 

2.49 

2.02 


624.43 

61.0 

-10 

89 

20 

1.31 

6.1 

73.9 

2.83 

0.30 

0.76 

8.8 

11.2 

10.8 

2.44 

2.44 

1.98 

(b) 

924.40 

61.0 

-17 

90 

55 

1.16 

6.1 

91.7 

2.53 

0.51 

0.84 

15.9 

20.5 

11.4 

6.77 

2.56 

2.02 


924.43 

61.0 

-17 

89 

20 

1.30 

6.1 

74.0 

2.83 

0.54 

0.76 

15.8 

20.4 

11.4 

2.46 

2.54 

2.00 


Figure 2. Effect of MVD on ice Shape. GLC 305 Airfoil. Ice-Shape Data from Tests by Chen 
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cm. Run 0- 
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(a) n = 0.4 
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Run 

c. 

ht, 

V, 

MVD, 

LWC, 

T, 

Po, 

A 

n 

u 

<P, 

e. 

Re, 

We, 

We c , 

We L , 

IvUIl 

cm 

°c 

m/s 

pm 

g/m 3 

min 

% 

A c 

u 

°C 

°c 

10 4 

10 3 

10 3 

10 6 

(a) 2-8-36 

91.4 

-12 

118 

44 

0.58 

14.6 

87.2 

2.67 

0.40 

0.57 

10.8 

11.8 

21.4 

9.54 

6.41 

5.28 

2-8-12 

30.5 

-13 

125 

26 

0.91 

3.0 

90.0 

2.68 

0.41 

0.55 

11.2 

11.8 

7.5 

6.29 

2.36 

1.96 

(b) 0-14-8-3 

91.4 

-15 

121 

38 

0.55 

13.5 

84.9 

2.40 

0.52 

0.54 

13.3 

14.8 

22.2 

8.59 

6.76 

5.54 

0-14-8-0 

30.5 

-15 

121 

20 

0.96 

2.6 

85.9 

2.42 

0.50 

0.54 

13.1 

14.5 

7.4 

4.54 

2.23 

1.83 


Figure 3. Effect of Model Size on ice Shape. GLC 305 Airfoil. Ice-Shape Data from Tests by Chen. 10 


decrease. Weber numbers with other lengths have been 
considered in the past, such as one based on droplet size 
and water properties, We: 

WeJ 2MVD ^ (14) 

o 

This Weber number seemed to be the logical choice 
based on the assumption that droplet splashing plays a 
role in establishing shapes for glaze ice. Anderson and 
Ruff 16 , Anderson 17 and Kind 3 have all used the Weber 
number of equation (14). Because in most studies, the 
MVD was scaled along with chord, matching We would 
lead to scale velocities not too different from those ob- 
tained from matching Wei, thus, positive scaling results 
from these studies were probably misleading. Weber 
numbers based on various water film thickness expres- 
sions have also been evaluated 18,19 with encouraging 
outcomes, and it is possible that L will prove to be a 
water-film thickness. 

With L °c d the scale velocity with We is = We LR is 



where the subscripts S and R refer to scale and refer- 
ence conditions, respectively. As a practical matter for 


tests with nearly constant air density, this result is little 
different from that obtained using constant We c (eq. 
(12)). The present study will provide preliminary results 
to suggest that equation (13) provides an effective simi- 
larity parameter to supplement Ruff s basic scaling 
method for SLD reference conditions. 


Test Description 


The scaling tests were performed in the NASA Glenn 
IRT. The IRT is a closed-loop, refrigerated, sea-level 
tunnel with a rectangular test section measuring 1.8 by 
2.7 m. It uses 10 spray bars to generate a cloud of su- 
percooled droplets. The Appendix-C cloud calibration 
used for these tests was performed in the summer of 
2000. 20 


The SLD calibration was made in the summer of 2002 
using the same methods as the Appendix C. Only a few 
specific MVD-LWC combinations for speeds of 51, 77 
and 103 m/s (100, 150 and 200 kt) have been calibrated 
to date. Therefore SLD tests are constrained to these 
particular conditions. 

The models used were NACA 0012 airfoil sections with 
chords of 91.4, 53.3 and 26.7 cm. The 91.4-cm-chord 
airfoil is pictured in figure 4 (a). It was a full-span, 
fiberglass model and served as the reference model. 
The 53.3-cm-chord model (1.7:1 scale) was of 61-cm 
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(a) 91.4-cm-Chord NACA 0012 Model Installed in IRT 
Test Section. 



(b) 53.3-cm-Chord NACA 0012 Model Installed in IRT 
Test Section. 

Figure 4. Model Description. 


span and made of aluminum. It was mounted vertically 
between splitter plates at the center of the IRT test sec- 
tion as shown in figure 4(b). Horizontal lines at the 
leading edge indicated tunnel center, ±2.5 cm and ±5 
cm from the center as guides for locating ice tracings. 
The 26.7-cm model (3.4:1 scale) was also of 61-cm 
span. Its appearance and mounting method were simi- 
lar to the 53.3-cm. All tests were run at 0° AOA al- 
though the mounting arrangement permitted rotation of 
the model for angle of attack changes. 

The IRT spray system has the ability to produce a stabi- 
lized spray within a few seconds of start; thus, shielding 
of the models was not needed during spray initiation. 

In preparing for a test, the temperature and airspeed in 
the test section and the air and water pressures on the 
spray manifolds were set. When these conditions had 
stabilized, the spray nozzle valves were opened to initi- 
ate the spray. The spray was timed for the required 
duration, and then turned off. The fan was brought to a 
full stop and the tunnel entered to record the ice shape. 
A heated ice knife with a cutout in the shape of the 
model was inserted into the leading edge of the ice to 
melt a thin slice down to the model surface. A card- 
board template was placed into this gap and an outline 
of the ice accretion traced. Tracings were taken at the 
vertical center of the tunnel (91 cm from the floor) and 
at 2.5 cm above the center. Ice shape differences be- 
tween the two tracing locations were never significant 
and only centerline shapes will be reported here. The 
tracings were digitized and the x-y coordinates for each 
ice shape recorded. The coordinates have been normal- 
ized by the model chord for presentation here. 

Uncertainty Analysis 

Uncertainties in the test conditions were estimated from 
differences between individual instalment readings, 
instrument uncertainty and fluctuations over the spray 
duration. Temperatures were believed to be good to 
±0.5°C, and the uncertainty in velocity was estimated to 
be 3%. For Appendix-C conditions the net uncertainty 
in MVD was estimated at ±12%. For SLD conditions it 
may have been as much as ±25%. These uncertainties 
are not referenced to an absolute value of MVD , which 
is unknown. Repeatability and scatter in the LWC cali- 
bration data suggests the uncertainty is about ±10% for 
both Appendix-C and SLD conditions. 

These uncertainties in the test parameters were esti- 
mated to produce the following uncertainties in the 
similarity parameters for the Appendix-C tests: 8% in 
! % , 10% in A c , 15% in n, 2% in Re, 13% in We and 5% 
in We L . For the SLD tests the uncertainties were: 4% 
in p 0 , 10% in A c , 15% in n, 2% in Re, 26% in We and 
5% in We,, 
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The effect of droplet size is shown in figure 5, which 
compares ice shapes obtained for SLD droplet sizes 
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MVD =40 pm, 10-12-02 Run 1 

(a) MVD, 200 and 40pm; c, 53.3 cm; V, 51 m/s; n, 0.3. 


A/I D =1 10pm. 10-12-02 Run 4 

MVD =20, 10-1 1-02 Run 3 

(b) MVD, 110 and 20pm; c, 53.3 cm; V, 51 m/s; n, 0.5. 


-0.04 0.00 0.04 0.08 0.12 0.16 

x/c 

MI D =95|im. 10-10-02 Run 2 
MVD =20|im. 10-10-02 Run 1 

(c) MVD, 95 and 20pm; c, 91.4 cm; V, 77 m/s; n, 0.3 


-0.04 0.00 0.04 0.08 0.12 0.16 


MVD =200pm. 10-10-02 Run 4 
MVD =20pm. 10-10-02 Run 1 

(d) MVD, 200 and 20pm; c, 91.4 cm; V, 77 m/s; n, 0.3 


Date/Run 


(b) 10-12-02/4 53.3 -11 51 


ht, 

v. 

MVD, 

LWC, 

°c 

m/s 

pm 

g/m 3 

-8 

51 

200 

1.20 

-6 

51 

40 

1.13 

-11 

51 

120 

1.03 


r, /3„, . , <p, 6, Re, We, We c , We L , 

min % c °C °C 10 4 10 3 10 3 10 6 

7.1 97.8 1.71 0.32 0.70 8.1 11.9 6.7 8.13 0.88 0.69 

10.8 84.0 2.44 0.28 0.57 6.1 9.0 6.6 1.62 0.87 0.68 

9.4 95.7 1.92 0.49 0.59 11.2 16.2 6.8 4.80 0.88 0.68 
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with those at Appendix-C conditions. Recorded test 
conditions and corresponding similarity parameters for 
each of the tests are given in a table accompanying the 
figure. Tests were run using the same NACA 0012 
model at the same velocity for each pair of shapes. 
Temperatures were determined to provide a match of 
the freezing fraction for each pair. The stagnation-line 
collection efficiency, /3 0 , could not be matched, but 
conditions were chosen so that the product of I3qA c 
would match. Later evaluation of calibration data indi- 
cated that the actual j3 f) A c for each pair of shapes proba- 
bly differed by 12 - 22%. This disagreement explains 
the differences observed in the leading-edge ice thick- 
ness for each pair of shapes. 

For figure 5(a), ice shapes obtained with 200 and 40pm 
are compared for the 53-cm model at 51 m/s. For 5(b) 
110 and 20pm are compared for the same model and 
velocity as for 5(a). Figures 5(c) and (d) show ice 
shapes recorded at 77 m/s for the 91 -cm model. 5(c) 
compares a 95-pm MVD shape with that of 20-pm, and 
(d) a 200-pm and 20-pm. For each of the comparisons, 
the main ice shape was little changed by the reduction 
in MVD from SLD conditions to an Appendix-C value. 

This correspondence of the main ice shape over wide 
ranges of drop size indicates that it should be relatively 
easy to simulate SLD conditions with Appendix-C 
droplet sizes, at least for the conditions studied here, if 
the main ice shape is of primary interest. However, two 
other features of those SLD accretions need to be noted. 
First, when drop sizes were larger than about 100pm 
the feather structures aft of the main ice shape were 
significantly larger than those of the Appendix-C 
shapes. These large feathers were particularly promi- 
nent for an MVD of 200pm. Further study is needed to 
determine under what conditions these growths occur, 
the physics behind them and whether they can be simu- 
lated in small-droplet accretions. 

Second, because the SLD collection efficiency was so 
much larger than that for Appendix-C conditions the 
SLD accretions featured small feathers well aft of the 
Appendix-C icing limits. Figure 6 gives recorded icing 
limits, normalized with respect to chord, for both the 
91- (solid symbols) and 53-cm models (open symbols) 
including all the Appendix-C and SLD test conditions 
for this study. These limits were recorded at the 
completion of each test at the same time the ice was 
traced. Each data point represents the average of 4 
estimates: the upper and lower surface limits for each 
of two ice tracing locations. Near the aft extent of 
icing, feathers become small and sparsely distributed. 
Therefore, significant variability in the estimated limit 
from run to run was experienced. The large scatter 
within the data for each model size was indicative of 



Po- % 

Figure 6. Icing Limits for Two NACA 0012 Models. 

for each model size was indicative of this variability. 
The results were estimated to have an uncertainty of 
about ±20%. 

The icing limits in figure 6 correlate with /3 n although 
the bulk of the data with the larger model tend to fall 
below those from the smaller model. A fit of the 53-cm 
data is shown with a solid line which passes through the 
origin. No correlation with test parameters was ob- 
served, indicating that it should be possible to deter- 
mine icing limits for SLD conditions from tests at Ap- 
pendix-C conditions if /3„ is matched. For scale tests it 
is always desirable to match //, to the reference value. 
However, if this cannot be done for the desired condi- 
tions, the main ice shape and icing limit can be deter- 
mined with separate tests: one with unmatched /3 n to 
find the main ice shape and one to determine the effect 
of /3 0 on icing limit. 

The results of figure 5 confirm Chen’s observations that 
main ice shapes are not affected by MVD. Thus, even 
for SLD conditions the We L of equation (13) is a more 
likely form of Weber number than the We of equation 
(14). To test the validity of finding scale velocity by 
matching scale and reference We L , SLD ice shapes were 
compared with those from Appendix-C conditions for 
which WeL matched. The results are given in figure 7. 
Reference ice shapes were obtained with a 91-cm-chord 
0012 model and an MVD of 160 pm. The velocity was 
51 m/s, LWC was 1.50 g/m 3 and the spray time 9.7 min. 
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Static temperature was varied from -11 to -19°C to give These freezing fractions were chosen to match values 

nominal freezing fractions of 0.3, 0.4 and 0.5. for previous Appendix-C tests using 53.3- and 26.7-cm- 



x/c 

c = 91.4 cm. 1-25-02 Run 1 
c = 53.3 cm. 1 1-13-00 Run 1 

(a) Model Chord Scaled from 91 to 53 cm; n, 0.3. 



x/c 

c- 91.4 cm. 1-25-02 Run 2 

c = 53.3 cm. 1 1-13-00 Run 4 

(b) Model Chord Scaled from 91 to 53 cm; n, 0.4. 



x/c 

c = 9 1 .4 cm. 2-8-02 Run 7 
c - 53.3 cm. 1 1-13-00 Run 5 

(c) Model Chord Scaled from 91 to 53 cm; n, 0.5. 



x/c 

c = 91.4 cm. 2-8-02 Run 7 
c = 26.7 cm. 1 1-16-00 Run 6 

(d) Model Chord Scaled from 91 to 27 cm; n, 0.5. 
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Figure 7. Ice-Shape Comparisons for Constant or Nearly Constant We^ 
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chord NACA 0012 models. 19 These Appendix-C ice 
shapes will serve as scale results for the present study. 
While it was not possible to match reference and scale 
Po, the values were within 12%, which as shown above 
is more than adequate for scaling main ice shapes. 
What is more important, the product PqA c matched for 
each pair of tests compared. Values of b, (p and 6 did 
not match. Conditions for the 53.3-cm model tests 
were such that We L matched that of the 91.4-cm, while 
We i for the 26.7-cm tests was about 30% higher than 
the reference for the example shown. The 26.7-cm- 
chord scale conditions providing an exact match of the 
reference We L will be tested in a future study. 


The generally good agreement of ice shapes in this 
study suggests that the parameters b, cp and 6, scale val- 
ues of which did not match the reference, can be ig- 
nored as long as the freezing fraction matches. 


Concluding Remarks 


Ice-shape evidence showed that an important similarity 
parameter needed in scaling analyses is the Weber 
number based on a length, L, proportional to model 
size. Thus, 



V-Lp v 

a 


(16) 


In figure 7 the shaded ice shape is that obtained with Good scaling was achieved in this study by matching 

the 9 1 .4-cm-chord at SLD conditions. The solid line scale and reference values of the parameters n and We,, 

shows the Appendix-C (scale) shape. Test conditions and the product /3oA c - Scale size ratios were as large as 

and corresponding similarity parameters for each set of 3.4:1 and freezing fractions covered the range from 0.3 

tests are given in the table accompanying the figure. to 0.5. 


Figure 7(a) shows the ice shapes for a freezing fraction 
of 0.3. The 53.3-cm-chord 38-pm-MFZ) test produced 
a main ice shape in very close agreement with the 91.4- 
cm-chord 160-qm-MFD shape. In addition, the sizes of 
both the large feathers adjacent to the main shape and 
the smaller feathers further aft were simulated well. 

Figure 7(b) compares the 91.4-cm reference with the 
53.3-cm-chord scale result for a freezing fraction of 0.4. 
The 53.3-cm-chord Appendix-C ice shape again gave 
an excellent match of the 91. 4-cm-chord SLD shape, 
although the lower-surface feathers just aft of the main 
shape were significantly larger for the SLD conditions. 

Figures 7(c) and (d) give results for a freezing fraction 
of 0.5 with the 91.4-cm reference conditions scaled to 
53.3 and 26.7 cm, respectively. Scale and reference ice 
shapes were in very good agreement for both scale 
models. The scale ice shapes for the 26.7-cm model 
matched the SLD reference for the freezing fractions of 
0.3 and 0.4, as well, showing that SLD ice shapes can 
be simulated at scale ratios of at least 3.4:1. 

Other Appendix-C tests with the 26.7-cm chord model 
were performed with different velocities than that 
shown in figure 7(d). None of these exactly matched 
the value satisfying constant We L \ however, the best 
match of the SLD shapes occurred when velocities 
close to the value required to match We, were selected. 

Appendix-C tests with these same conditions have been 
repeated in several other tunnel entries over the last 2 
years. Because ice shapes vary somewhat from one 
entry to another, agreement of Appendix-C and SLD 
results for other entries was not always as good as 
shown here. 


Good agreement of both quantity and shape of the main 
ice accumulation between Appendix-C scale conditions 
and 1 60- um SLD reference conditions was observed. 
Tests that directly compared ice shapes formed at 
200(im and at 20pm also showed good agreement of 
the main accretions. These results argue against a sig- 
nificant role of splashing for the SLD conditions con- 
sidered. However, the mechanism of formation of large 
feather structures seen in SLD accretions for MVD’s 
larger than about 100pm needs to be researched. The 
present SLD reference tests were made with velocities 
of about 50 and 77 m/s, and these conclusions may not 
be valid for higher velocities. 

Additional analysis is needed to identify the appropriate 
length L for use in the Weber number of equation (16). 
The ice shapes analyzed and compared in this study 
suggest that L is either independent of MVD and L WC, 
or only weakly dependent, while it is proportional to 
model size. 

This and other scaling studies have shown that LWC 
can be chosen fairly arbitrarily as long as reference and 
scale freezing fractions match. Excellent scaling results 
were obtained in this study without matching b, tp, or 6. 
It should be possible to gain flexibility in scaling tests, 
then, by choosing convenient values for LWC and tem- 
perature with the only constraint being that the freezing 
fraction matches the reference value. It was also shown 
that, while it is desirable to match fio, main ice shapes 
appear to be independent of MVD. The ability to set 
scale MVD to values convenient to the test facility adds 
additional flexibility in determining scale conditions. 
The main effect of MVD (through p n ) is on the icing 
limit, but because icing limits can be correlated with /fy 
tests of scale-model icing limit for a range of fio can be 
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run independently of those recording the main ice 
shape. This approach should provide adequate informa- 
tion about that characteristic. 
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